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ABSTRACT 


I 

\ 

This  paper  is  concerned  with  piecewise  smooth  plane 
deformations  in  an  isotropic,  incompressible  elastic  material. 
An  explicit  necessary  and  sufficient  condition  for  the  existence 
of  piecewise  homogeneous  equilibrium  states  is  established,  and 
the  set  of  all  such  states  is  precisely  characterized.  A 
particularly  simple  expression  is  derived  for  the  ^driving 
traction'^  on  a  surface  of  discontinuity  in  the  deformation 
gradient. 
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1 .  Introduction 

When  finitely  deformed,  certain  elastic  materials  are  capable  of 
sustaining  deformation  fields  which  involve  discontinuous  gradients.  In 
such  circumstances,  the  deformation  varies  smoothly  throughout  the  body 
except  at  certain  surfaces  across  which  the  displacement  gradient  suffers  a 
jump  discontinuity,  even  though  the  displacement  itself  remains  continuous. 
Continuum  mechanical  treatments  of  stress  - induced  phase  transformations  in 
solids  involve  such  deformations;  see,  for  example,  [1-3]. 

In  this  paper  we  first  examine  conditions  under  which  a  piecewise 
homogeneous  plane  deformation  can  be  sustained  by  an  incompressible,  iso¬ 
tropic  elastic  material.  The  question  of  the  existence  of  such  fields  has 
been  investigated  previously.  For  example,  the  analysis  in  [4]  addresses 
this  issue  within  the  context  of  compressible  elastic  materials  undergoing 
plane  deformations,  while  [5]  contains  a  similar  study  for  incompressible 
materials,  and  the  investigation  in  [6]  approaches  the  subject  from  a  vari¬ 
ational  point  of  view.  By  restricting  attention  in  this  study  to  materials 
that  are  isotropic  and  incompressible  and  to  plane  deformations,  we  are 
able  to  derive  a  single  necessary  and  sufficient  condition  for  the 
existence  of  piecewise  homogeneous  deformations.  Moreover,  this  condition 
is  expressed  in  a  form  that  is  particularly  useful;  in  addition  to  provid¬ 
ing  information  on  existence,  it  also  allows  us  to  characterize  the  set  of 
all  possible  piecewise  homogeneous  equilibrium  states.  A  similar  result 
for  anti -plane  deformations  was  established  in  [7], 

When  the  theory  of  elasticity  is  broadened  to  allow  for  equilibrium 
fields  with  discontinuous  deformation  gradients,  it  turns  out  that  the 
usual  balance  between  the  rate  of  external  work  and  the  rate  of  storage  of 
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elastic  energy  during  a  quasi-static  motion  no  longer  holds;  see  [8] 
Instead,  one  finds  that  mechanical  energy  may  be  dissipated  at  points  on 
the  surfaces  of  discontinuity.  This  is  analogous  to  the  dissipation  of 
energy  in  an  invisc id  compressible  fluid  when  the  flow  involves  a  shock; 
the  analogy  suggests  the  term  "equilibrium  shock"  for  surfaces  bearing 
discontinuities  in  the  displacement  gradient  in  elastostatic  fields. 

As  discussed  in  [9-12],  the  altered  energetics  of  fields  with  equili¬ 
brium  shocks  lead  to  the  notion  of  a  scalar  "shock  driving  traction"  f 
which  may  be  viewed  as  a  normal  traction  that  the  body  applies  to  the  shock 
at  each  of  its  points.  At  each  point  on  a  moving  shock  during  a  quasi - 
static  motion,  the  product  of  f  with  the  normal  velocity  of  the  shock  rep¬ 
resents  the  local  rate  of  dissipation.  The  presence  of  dissipation  suggests 
that  elastic  fields  with  equilibrium  shocks  might  be  used  to  model  certain 
types  of  inelastic  behavior  in  solids.  Investigations  that  confirm  this 
possibility  are  summarized  in  [10-12],  where  the  nature  of  the  shock  trac¬ 
tion  f  plays  a  major  role. 


The  problems  treated  in  [10-12]  are  essentially  one  -  dimens ional ;  for 
these,  the  structure  of  the  shock  traction  and  its  dependence  on  the  defor¬ 
mation  are  well  understood.  For  problems  involving  more  dimensions,  the 
matter  is  less  fully  explored,  although  Yatomi  and  Nishimura  [13]  have 
obtained  some  valuable  results.  The  second  principal  objective  of  the  pre¬ 
sent  paper  is  the  derivation  of  an  especially  simple  expression  for  the 
shock  driving  traction  for  isotropic,  incompressible  elastic  materials  in 
the  case  of  plane  deformations. 


In  Section  2  we  present  some  background  material  on  equilibrium  shocks 


in  the  appropriate  setting.  The  structure  of  piecewise  homogeneous 
equilibrium  states  is  described  in  Section  3,  and  a  necessary  and 
sufficient  condition  for  the  existence  of  such  states  is  established  and 
discussed  in  Section  4.  In  Section  5,  the  representation  for  shock 
traction  mentioned  above  is  derived  and  discussed. 


2 .  Preliminaries  on  Finite  Plane  Strain  and  Equilibrium  Shocks. 

Consider  an  incompressible  body  which,  in  a  reference  configuration, 
occupies  a  cylindrical  region  of  space;  denote  by  D  the  open  middle 
cross-section  of  this  region.  Let  X  -  (  0;  ,  £2  >  S3 1  be  a  fi*ed 

rectangular  cartesian  frame  whose  origin  0  and  basis  vectors  e i ,  £2  are 
the  plane  of  D. 

A  plane  deformation  of  the  body  is  described  by  the  one-to-one  mapping 
y  ”  y(x)  -  x  +  u(x) ,  where  x  and  y  are  the  position  vectors  of  a  particle 

A 

in  D  and  D*  -  y(D) ,  respectively.  For  the  present,  it  is  assumed  that 
u  e  C^(D),  and  we  let  F  and  £  denote  the  two-dimensional  deformation 
gradient  tensor  and  the  left  Cauchy-Green  deformation  tensor,  respectively: 


F  -  Vy,  G  -  FFT . 

(2. 

1  ) 

Since 

the  material  is  incompressible, 

the  deformation  invariants  I 

and 

J 

obey 

1  “  tr(fXT)  ”  Al2  +  *22  5  2’ 

J  -  det 

F  -  XjA2  -  1, 

(2 

2) 

where 

A^(x^)  >  0  and  ^2^10  -  0 

.  are  the 

principal  stretches 

of 

the 

deformation  at  a  point  in  in  D. 


..5-- 


The  following  kinematical  result  has  been  established  in  [5]:  Every 


two-dimensional  tensor  £  with  unit  determinant  admits  the  two  decomposi¬ 


tions 


l  -  2ixi  “  £222. 


(2.3) 


where  >  Q2  are  ProPer  orthogonal  tensors,  and  K]_ ,  £2  >  are  tensors  whose 


matrices  of  components  in  suitable  rectangular  cartesian  frames  and  X2 


have  the  common  value 


[KXXI ]  -  [K2X*] 


(2.4) 


k  -  7(1-2)  >  0,  where  I  -  trFFj. 


(2.5) 


It  follows  that  a 


lane  volume -preserving  deformation  can  be  decomposed 


simple  shea 


a  _ suitable  direction 


rotation:  k  denotes  the  amount  of  this 


shear . 


Next,  let  o(x)  be  the  two-dimensional  nominal  stress  tensor  field  on  D 


and  the  corresponding  Cauchy  stress  tensor  field  on  D* : 


o(x)  -  r( y(x))  F'T(x) ,  xeD . 


(2.6) 


If  C  is  a  smooth  arc  in  D,  and  C*-^(C)  is  its  image  after  deformation,  the 


nominal  traction  ^  on  C  and  the  true  traction  jt  on  C*  are 


s  -  crN  on  C , 


t  -  rn  on  C* , 


WS 
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where  N  and  n  are  corresponding  unit  normals  to  C  and  C*  respectively. 
Assuming  for  the  present  that  o(x)  is  continuously  differentiable  on  D, 
equilibrium  in  the  absence  of  body  forces  requires 


div  £7-0, 


or  equivalently, 


div  r  -  0 , 


o FT  -  FctT 


on  D, 


(2.8) 


on  D* . 


(2.9) 


Turning  to  the  constitutive  law  of  the  incompressible  material  at  hand, 
we  suppose  that  it  is  homogeneous,  isotropic  and  hyperelastic.  The  plane 

A 

strain  elastic  potential  W(F)  then  has  the  form 


W(F)  -  W(I) ; 


(2.10) 


the  value  of  W  is  the  strain  energy  per  unit  reference  volume.  The  associ¬ 
ated  constitutive  law  takes  the  equivalent  alternate  forms 


a  -  2W'(I)F  -  pF*  ,  r  -  2W'(I)FFT  -  pi, 


(2.11) 


where  the  scalar  field  p  arises  because  of  the  incompressibility  con¬ 
straint  . 


The  particular  deformation  described  by 


yi  "  x!  +  kx2- 


y2  -  x2, 


(2.12) 


is  a  simple  shear,  where  the  constant  k  denotes  the  amount  of  shear.  The 
corresponding  shear  stress  component  r^2  found  from  (2.12),  (2.11), 
(2.1),  (2,2)  to  be  r^2  “  r*(k)  ,  where 


r .  /,  /,  W*  A/,  f-  ,  .  /.r  X,  f , *.  <*.  / .  •  .  ..  .  *  .  •  ,  ik.  • 


t (k)  -  2kW' (2+kz) , 


<  k  <  ® . 


(2.13) 


The  odd  function  r(k)  is  the  shear  stress  response  function  of  the  material 
in  simple  shear.  Throughout  thi  paper  it  will  be  assumed  that  the  secant 
modulus  of  shear  M(k)  is  positive: 


M(k)  -  2  W'(2+k2)  >0  for  -®  <  k  <  ®  ; 


(2.14) 


this  is  consistent  with  the  Baker- Ericksen  inequality.  While  materials  for 
which  M(k)  is  not  always  positive  are  of  interest  (see  [11]),  this  case 
will  not  be  pursued  here.  Observe  from  (2.13)  that 


7(1-2) 


r (k)  dk  for  I  >  2, 


(2.15) 


which  implies  that  the  in-plane  response,  in  all  plane  deformations,  of  an 


isotropic 


essible  elastic  material  is  completely  characterized  by 


The  displacement  equations  of  equilibrium  for  plane  strain  in  a  homoge¬ 
neous,  isotropic,  incompressible  material  are,  by  (2.11),  (2.8),  (2.2), 


capy8  (£) uy  ,  fiS  "  P,/3  ® 


det  F 


on  D , 


(2.16) 


where  FQ/3  -  fiaj3  +  ua  /3  and 


ca0y6 (£>  “  2W'(I)6Q75^{  +  4W"  (I)Fa^F7fi . 


(2.17) 


Here,  a  comma  followed  by  a  subscript  indicates  partial  differentiation 
with  respect  to  the  corresponding  x- coordinate ,  while  repeated  subscripts 


are  summed  over  1,2. 


It  has  been  shown  in  [5]  that,  in  the  presence  of  (2.14),  the  system  of 
partial  differential  equations  (2.16)  is  elliptic  at  a  solution  u,  p  and  at 
a  point  x,  if  and  only  if 

r ' (k(x))  >  0,  (2.18) 

where  k(x)  -  7(I(x)-2)  -  7(trJF(x) F^(x) -2)  is  the  amount  of  shear  associated 
with  the  deformation  at  that  point.  From  (2.13),  (2.14)  one  has  t’(0)  >  0, 
so  that  ellipticity  necessarily  prevails  at  infinitesimal  deformations.  If, 
however,  r  fails  to  be  monotonically  increasing  on  -®  <  k  <  °°,  ellipticity 
will  be  lost  if  the  deformation  is  severe  enough.  If  r ' (k)  >  0  for  all  k, 
we  say  that  the  material  is  elliptic  for  plane  deformations. 

One  consequence  of  a  loss  of  ellipticity  of  the  governing  partial  dif¬ 
ferential  equations  is  the  possible  occurrence  of  elastostatic  fields  which 
are  less  than  classically  smooth.  In  order  to  account  for  such  weak  solu¬ 
tions,  the  smoothness  assumptions  made  previously  must  be  relaxed.  Of  par¬ 
ticular  interest  is  the  case  wherein  the  field  quantities  possess  the 
aforementioned  degree  of  smoothness  everywhere  expect  on  one  or  more  arcs 
in  D.  Accordingly,  we  now  allow  for  the  possibility  that  although  u  is  con¬ 
tinuous  in  D,  there  is  a  single  smooth  curve  C  c  D  such  that  a  and  u  are 
respectively  once  and  twice  continuously  differentiable  in  D-C,  and  £  and 
Vu  suffer  finite  jump  discontinuities  across  C.  In  these  new  circum¬ 
stances,  the  field  equations  discussed  previously  continue  to  hold  in  D-C. 
In  addition,  equilibrium  considerations  require  that  the  nominal  traction  s 
be  continuous  across  C;  this  in  turn  implies  that  the  true  traction  t^  is 
continuous  across  C*,  the  deformation- image  of  C.  Thus 


[ a  1 ] N  -  0  on  C , 


I  t  1 ] n  -  0  on  C* , 


(2.19) 


where  N  and  n  are  unit  normals  on  C  and  C*  respectively,  and  [[•]]  indi¬ 
cates  the  jump  across  the  appropriate  curve.  Note  that  displacement  conti¬ 
nuity  across  C  requires 


[  F] ]L  -  0  on  C, 


(2.20) 


where  L  is  a  unit  tangent  vector  on  C .  A.  curve  C  carrying  jump  discontinu¬ 
ities  in  F,  p  and  j?  while  preserving  continuity  of  displacement  and  trac¬ 
tion  is  called  an  equilibrium  shock. 


3 .  Piecewise  homogeneous  equilibrium  states. 

In  order  to  investigate  many  of  the  local  issues  related  to  equilibrium 
shocks,  it  is  sufficient  to  consider  the  case  in  which  D  coincides  with  the 
entire  (x^ , X2) -plane ,  C  is  a  straight  line  through  the  origin,  and  the 
deformation  gradient  £  as  well  as  the  hydrostatic  pressure  p  are  constant 
on  either  side  of  C.  Under  these  conditions,  the  field  equations  (2.8), 
(2.11)  will  be  trivially  satisfied  in  D  -  C,  and  the  requirements  remaining 
to  be  fulfilled  are  (2.20)  and  either  of  (2.19). 


Let  L  be  a  unit  vector  along  the  straight  line  C  through  the  origin, 

and  let  N  —  ei  x  L  be  the  unit  normal  to  C  associated  with  L,  so  that 

-+  . 

(L.N.jirj)  is  a  right-handed  orthonormal  triplet.  Let  D  -  (x  |  x*N  >  0), 

D  -  |x  |  x*N  <  0)  be  the  two  open  half-planes  into  which  C  divides  D;  see 

■f* 

Figure  1.  Points  in  D  lie  on  the  pos i t i ve  side  of  the  shock 


Consider  the  piecewise  homogeneous  deformation 


F  x  for  xeD  , 


F  x  for  xeD  , 


where  F  and  F  are  distinct,  constant  unimodular  tensors 


(3.1) 


F  *  F  , 


det  F  -  det  F  *  1 


(3.2) 


+  —  +  “ 

Define  G,  G,  I,  I,  k.  and  k  by 

+  +  + 

G  -  F  FT  , 


tr  G  , 


k+  -  (I  -  2)1/2  >  0  ,  k_  -  (I  -  2)1/2  >  0 


(3.3) 


(3.4) 


(3.5) 


In  view  of  the  kinematical  result  (2.  3) -(2. 5),  k+  and  k.  represent  th 
respective  amounts  of  shear  associated  with  the  deformations  on  D  and  D  . 


The  displacement  corresponding  to  the  deformation  (3.1)  will  be  conti 
nuous  across  C  if  and  only  if 


F  L  -  F  L 


Assume  that  (3.6)  holds,  and  let 


(3.6) 


l  -  g  FL  -  g  FL  . 


(3.7) 


where 


g  -  1/| FL|  -  1/ I FL |  >  0; 


(3.8) 


I  1*1  IS  IV 


--li¬ 


the  unit  vector  is  tangent  to  the  image  C*  of  C  under  the  deformation 


(3.1).  The  vector  n  defined  by 


n  -  e-3  x  £ 


(3.9) 


is  the  unit  normal  to  C*  associated  with  (£,n,<|3)  is  a  right-handed 

orthonormal  triplet.  Using  the  definitions  of  N,  2  and  n  as  well  as  the 

facts  that  det  F  >  0,  det  F  >  0,  one  can  show  that  D*  —  (y  |  y • n  >  0)  and 

D*  -  (y  |  y*n  <  0)  are  the  respective  images  of  D  and  D  under  the 

deformation  (3.1).  Because  det  F  -  det  F  -  1,  (3.7)  and  (3.9)  imply  that 


FTn  -  FTn  =  g  N, 


from  which  it  further  follows  that 


(3.10) 


g  -  (n  •  Gn)!/^  _  (n  .  Gn) 1/2 . 


(3.11) 


Note  from  (3.11)  that 


[[n  •  Gn]]  -  0, 


(3.12) 


where  [[•]]  indicates  the  jump  across  the  shock. 


A  fiber  whose  orientation  in  the  reference  configuration  is  parallel  to 
the  shock  C  suffers  a  relative  extension  (l-g)/g  ;  when  g  -  1,  such  a  fiber 
is  unextended  by  the  deformation  (3.1).  When  this  is  the  case,  we  say  the 
shock  is  normal :  otherwise,  it  is  oblique . 


Let  a  unit  vector  L  and  a  unimodular  tensor  ^F  be  given.  Define  g  and  2 
through  the  first  equalities  in  (3.8)  and  (3.7),  respectively,  and  define^ 
by  (3.9).  A  unimodular  tensor  F  will  satisfy  (3.6)  if  and  only  if  there  is 
a  number  k  such  that 


-  - 12  -  - 


+ 

F  -  (1  +  *  t  ®  n)F  ,  (3.13) 

/sj  <V  'X-  ^ 

where  ®  denotes  tensor  multiplication  of  two  vectors.  To  prove  this,  sup- 

—  "t 

pose  first  that  F  satisfies  (3.6).  It  follows  that  £  -  F  +  h  ®  N  ,  where  N 

+  -  -r  i  + 

-  y  x  L  and  h  -  FN  -  FN.  Alternatively,  we  may  write  F  -  (1_  +  (heNJF'1)?  - 

+  i  T  +  - 

(1  +  h®(F" L) 1N) F  so  that  by  the  first  equality  in  (3.10),  F  =» 

+ 

(1_  +  (l/g)h  ®n)F.  Taking  the  determinant  of  both  sides  of  this  equality  and 
recalling  (3.2)2  implies  that  h  •  n  -  0.  Consequently  h_  is  parallel  to  2, 
and  we  may  write  h  —  (*g)  £*,  which  establishes  (3.13)  .  The  definition  of  h 
together  with  (3.10)  and  k  -  (l/g)h*2  gives 

k  -  -  (1/g2) [ [i_  •  Gn] ]  .  (3.14) 

Conversely,  if  F  is  given  by  (3.13),  tnen  it  is  unimodular  and 
satisfies  (3.6)  . 


V/  v  vdv:-.-:  ^ 


AN 
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while  the  dot  product  of  (3.15)  with  l  gives 


l  •  [ [M(k)G] ]  n  -  0  . 


(3.17) 


Conversely,  if  (3.16)  and  (3.17)  hold,  then  so  does  (3.15).  Finally,  using 


(3.14)  to  eliminate  i»Gn  in  (3.17)  yields 


*  “  [M(k+)  -  M(k_)]  (i,  •  Gn)/  [g2M(k.) 


(3.18) 


Several  useful  facts  follow  from  (3.18).  First,  in  a  piecewise  homogeneous 


equilibrium  state,  necessarily  k+  *  k_ ;  otherwise,  by  (3.18),  k  -  0,  so 


that  by  (3.13),  F  -  F  ,  contradicting  the  first  of  (3.2).  Second,  one  must 


have  k+  +  0,  k_  *  0.  If,  for  example,  k+  -  0,  then  (3.3)  and  (2.3),  (2.4) 

+•  + 

applied  to  F  show  that  G  is  the  identity.  From  (3.18)  it  again  follows 


that  k  —  0,  which  is  impossible.  Similarly,  k.  r*  0.  Thus  the  deformations 


on  both  sides  of  the  shock  must  involve  nontrivial,  distinct  amounts  of 


shear  k+,  k.  . 


Let  a  unimodular  tensor  F  be  given.  The  shock  problem  requires  the 


determination  of  a  unimodular  tensor  F  +  F  and  a  unit  vector  L  such  that 


(3.6)  and  (3.17)  hold;  in  (3.17),  G  and  G  are  to  be  given  by  (3.3),  k+  and 


,+  ,  +  + 

k.  by  (3.5),  g  -  1/|FL| ,  i  -  g  FL  ,  and  n  is  as  in  (3.9).  For  a  given  F, 


let  F,  L  be  a  solution  of  the  shock  problem,  and  let  p  be  a  given  constant. 


Define  p  by  means  of  (3.16);  then  F,  F\  p,  and  p  will  generate  a  piecewise 


homogeneous  equilibrium  state  through  (3.1)  and  the  constitutive  law 


(2.11)]^.  The  shock  C  is  the  straight  line  through  the  origin  determined  by 


m 
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4 .  Solutions  of  the  Shock  Problem. 

Here  we  establish  a  necessary  and  sufficient  condition  for  the 
existence  of  solutions  to  the  shock  problem,  and  we  discuss  its  interpreta¬ 
tion  and  some  of  its  implications . 


Proposition:  Let  F  be  a  given  unimodular  tensor,  and  let  k+  be  given 
by 


k+  -  7( tr(F  Ft)  -  2}  . 


(4.1) 


A  solution  to  the  shock  problem  corresponding  to  F  exists  if  and 
only  if  there  is  a  number  k.  *  k+  such  that 


[r(k+)  -  r(k_)](k+  -  k.)  <  0  , 


(4.2) 


where  r(k)  is  the  shear  stress  response  function  characteristic  of 
the  material. 


We  first  assume  that  the  shock  problem  has  a  solution  F,  L  for  the 
+• 

given  F  and  show  that  the  result  stated  above  is  necessary.  We  begin  by 
recording  the  identity 


(n*Gn)2  +  (i-Gn)2  -  |Gn|2  -  (n*G2n) . 

IV  <V/V  1  /v' 


(4.3) 


The  Cayley -Hamilton  theorem  applied  to  the  unimodular  tensors  G  and  G, 
together  with  (4.3)  and  (3.11),  yields 


+ 

r»4  ,  1-2 


(!-Gn)2  ~  -g4  +  Ig2  -  1  ,  (£'Gn)2  -  -g4  +  Ig2  -  1  . 


(4.4) 


Using  (4.4)  and  (3.5)  in  (3.17)  gives  the  following  equation  for  g  in  terms 
of  k+  and  k.  : 


[[M2(k)]]+g4  *  [[M2(k)(2+k2)]]+g2  +  [ [M2 (k) ] ] +  -  0 


(4.5) 


Because  of  ( 3 . 11)^,  we  may  write 


.'.vV a.ww'v'V, V-V-'.--'.--' 


TiV 


»  ►  -S  W  V  j 
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g 


2 


+  + 

X^2  cos20  +  >2^  sin2#  , 


(4.6) 


where  9  is  the  angle  in  [0,w]  between  n  and  the  principal  direction  of 

that  corresponds  to  the  principal  stretch  X^.  Without  loss  of  generality, 

■+■  +- 

we  assume  that  X^  <  X2 ,  so  that  it  then  follows  from  (4.6)  that 


+  + 


Xi  <  g  <  x2. 


(4.7) 


0  "*■ 

Returning  to  (4.5),  we  note  that  if  [[Mz(k)]]_-  0,  we  would  have  g  -  0, 
which  is  impossible.  Thus  (4.5)  may  be  rewritten  as 

(fM2(k)k2]]+/  [ [M2 (k) ] ] +  -  g2  +  g'2  -  2  >  0.  (4.8) 

Since  M(k)  -  r(k)/k  ,  (4.8)  is  equivalent  to 

[ [ t2(k) ] ]  /  [ [ r2 (k)/k2 ] ]  -  g2  +  g'2  -  2  >  0.  (4.9) 

+  +  o  T 

From  (4.6)  and  the  fact  that  X^  -  I/X2 ,  one  shows  readily  that  gz  +  g'z  < 
■f  9  +  9  +  9 

Xf  +  X2  -  I  —  2  +  k+z .  It  then  follows  from  (4.9)  that  necessarily 

0  <  [[r2(k)]]^  /  [ [r2(k)/k2] ]+<  k+2  .  (4.10) 

Finally,  one  can  show  that  (4.10)  is  equivalent  to 

[r(k+)  -  r(k_)](k+  -  k.)  <  0  ,  (4.11) 

thus  establishing  the  necessity  of  (4.2). 

We  now  show  that  the  conditions  stated  in  the  proposition  are  suffi- 
cient  for  the  existence  of  an  equilibrium  shock.  Let  F  be  a  given  unimodu- 
lar  tensor,  and  let  k+  be  defined  by  (4.1).  Suppose  that  there  is  a  number 
k.  *  k+  satisfying  (4.2);  note  that  there  may  be  many  such  numbers.  There 


fo-+ 


is  then  a  root  g  of  (4.9)  satisfying  (4.7)  for  this  k.  ;  in  fact,  there  are 
two  such  roots.  For  each  such  g,  there  is  an  angle  6  satisfying  (4.6)  and 

hence  determining  a  unit  vector  n_.  Corresponding  to  this  pair  g,n_,  one 

•*-T 

defines  N  -  (1/g)  FAn;  see  (3.10).  Let  L  *  N  x  ej,  £  “  E  x  ®3 ■  and  define  k 

+■  +  Tt  -  - 

by  (3.18),  using  G  -  F  F  .  Finally,  define  F  by  (3.13).  With  L  and  con¬ 
structed  in  this  way,  (3.6)  is  automatically  satisfied.  It  is  readily  veri¬ 
fied  that  (3.17)  holds  as  well.  Thus  F  and  L  furnish  a  solution  of  the 

+ 

shock  problem  corresponding  to  the  given  F,  and  the  proof  is  complete. 

The  proposition  just  proved  may  be  interpreted  geometrically  with  the 
help  of  a  set  £  in  the  (k^.k2) -plane  defined  as  follows: 

£  -  { (kx ,k2)  |  [r(ki)  -  r(k2)](k1  -  k2)  <  0,  kx>0 ,  k2>0 ,  kL  *  k2 } .  (4.12) 

+ 

Given  a  unimodular  tensor  F,  the  proposition  asserts  that  the  corresponding 
shock  problem  has  a  solution  if  and  only  if  there  is  a  number  k_  such  that 

(k+,k_)  e  E,  where  k+  -  y[tr(FFA)  -  2]  is  the  amount  of  shear  associated 

+ 

with  F.  For  illustration,  consider  a  material  whose  shear  stress  response 

A 

function  r(k)  is  as  shown  in  Figure  2.  Let  T  be  the  closed  curve  defined  by 

r  -  l(kltk2)  I  r(kx)  -  r(k2),  kj_  >  0,  k2  >  01;  (4.13) 

the  set  E  for  this  material  consists  of  all  points  on  and  within  F  except 
for  the  points  on  the  diagonal  line  k^  -  k2 .  E  is  symmetric  about  this 
diagonal.  Figure  3  shows  the  curve  T  and  the  region  E  for  a  material  of  the 
kind  described  above  whose  shear  stress  response  is  given  by 

r(k)  -  ji(k  -  yk^  +  k^)  ,  7  -  1.875  ,  n  >  0  ,  (4.14) 

the  qualitative  nature  of  Figure  3,  however,  does  not  depend  on  this  spe¬ 
cial  choice,  as  long  as  r(k)  has  the  form  shown  in  Figure  2 


The  shaded 


A  A 

region  in  Figure  3  consists  of  points  (k^,k2>  at  which  r'(k^)  and  r'(k2; 
are  both  positive;  the  corresponding  piecewise  homogeneous  equilibri  in 
states  always  involve  deformations  that  are  elliptic  on  both  sides  of  the 
shock.  This  is  not  the  case  for  points  (k^,k2>  in  the  unshaded  part  of  Z. 

For  the  special  case  of  a  normal  shock,  g  -  1,  and  (A. 9)  immediately 
yields 

r'Ck+l  -  t  (k.)  ;  (4.15) 

thus  the  set  of  points  in  E  that  correspond  to  normal  shocks  are  precisely 
the  points  on  the  boundary  T  of  E  except,  of  course,  for  the  points  with  k  + 
-  k.  . 

A 

If  r’(k)  >  0  for  all  k,  so  that  the  shear  stress  in  simple  shear  is  a 

monotone  increasing  function  of  the  amount  of  shear,  the  material  under 
consideration  is  elliptic  for  plane  deformations.  It  is  an  immediate 

consequence  of  the  proposition  that  E  is  empty  in  this  case,  and  no  piece - 

wise  homogeneous  equilibrium  states  of  plane  strain  exist  for  such  a  mate¬ 
rial.  On  the  other  hand,  if  the  material  is  such  that  r ' (k)  is  non-positive 
throughout  some  interval,  then  -  since  r  ’  (0)  >  0  by  (2  13),  (2.14)  -  piece - 

wise  homogeneous  plane  equilibrium  deformations  will  certainly  exist  for 

+• 

suitably  chosen  values  of  F. 

The  observations  of  the  preceding  paragraph  were  made  in  [5'  on  the 
basis  of  a  different  analysis.  In  addition  to  providing  information  on  the 
existence  of  shocks,  the  approach  taken  here  also  characterizes  all  pos¬ 
sible  piecewise  homogeneous  equilibrium  states  corresponding  to  a  given  F. 
For  each  k  >  0,  let  E^  be  the  set  in  the  (kj , kj )• plane  defined  by 


-  {*  |  (k,/t)  6  SI  ■  u  |  (x,k)  e  Z) 


(4.16) 


Zjj  is  the  cross-section  of  Z  at  either  kj  -  k  or  k2  -  k  .  If  £  is  a 
unimodular  tensor,  let  k  -  J[ tr(FFj)  -  2];  if  Z^  is  non-empty,  then  there 
exists  an  equilibrium  shock  corresponding  to  every  k.  in  Z^  . 

We  turn  now  to  a  discussion  of  the  structure  of  the  driving  traction 
associated  with  equilibrium  shocks. 


5.  Shock  Traction 

For  the  purposes  of  this  section,  it  is  necessary  to  consider  quasi - 
static  motions  of  the  body.  Let  (u(«,t),  p(*,t)|,  t0  <  t  <  t^ ,  be  a 

one -parameter  family  of  weak  solutions  of  the  displacement  equations  of 
equilibrium  (2.16)  of  the  type  described  in  Section  2;  let  Ct  c  D  be  the 
associated  family  of  shocks.  Assume  that  the  particle  velocity  v(x,t) 
du(x,t)/dt  exists  and  is  continuous  in  (x,t)  for  x  e  D-Ct,  tg  <  t  <  t^ , 
and  that  v  is  piecewise  continuous  on  D  x  [tg.tjJ. 


Let  V(x,t)  be  the  velocity  of  a  point  on  the  moving  curve  Ct  which 
at  time  t  is  located  at  x,  and  denote  by  N(x,t)  the  unit  normal  to  Ct  at 
this  same  point.  Note  that,  even  though  the  deformations  here  need  not  be 
piecewise  homogeneous,  equations  (3.2)-(3.18)  and  (4.1)-(4.12)  continue 

to  hold  at  each  point  on  Ct  and  each  instant  t,  provided  that  one 

♦ 

interprets  F(x,t)  or  F(x,t)  as  the  limiting  values  of  F  as  the  point  x  on 
Ct  is  approached  from  its  positive  or  negative  side,  respectively.  The 
positive  side  of  Ct  is  the  one  into  which  ^  points . 


mmmm 
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As  shown  in  [8],  the  presence  of  the  shock  affects  the  balance  of 
mechanical  energy  during  the  motion  in  an  essential  way.  For  any  fixed 
regular  sub-region  n  of  D,  let  d(t)  represent  the  difference  between  the 
rate  of  external  work  on  IT  and  the  rate  at  which  elastic  energy  is  being 
stored: 


d(t) 


Js*v  ds 

an 


d_ 

dt 


W(F)  dA, 

n 


tQ  <  t  <  ti: 


(5.1) 


d(t)  is  the  rate  of  dissipation  of  mechanical  energy  associated  with  the 
particles  which  occupy  the  region  n  in  the  reference  state.  By  adapting  to 
plane  deformations  an  argument  given  in  [8],  one  can  show  that  d(t)  may  be 
written  as 


d(t)  -  | 


[(P])*n  -v  ds, 


where  P(x,t)  is  the  two-dimensional  energy -momentum  tensor, 


(5.2) 


l  "  w(£)i  -  £t  °>  x  e  D-ct,  t0  <  t  <  tL.  (5.3) 

If  £(»,t)  happens  to  be  continuous  across  Ct  at  an  instant  t,  (5.2) 
gives  d(t)  -  0.  In  general,  however,  if  II  intersects  Ct,  d(t)  *  0  . 


It  has  been  shown  in  [8]  that,  because  of  displacement  and  traction 
continuity  across  Ct  ,  the  vector  [ [P] ]  N  is  normal  to  Ct  at  each  of  its 
points.  Thus  we  may  write  (5.2)  in  the  alternate  form 


d(t) 


f  N*V  ds, 


Ctnn 


where  f(x,t)  is  defined  by 


(5. A) 
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f  -[[£]]  JS  •  N,  X  e  ct,  t0  <  t  <  t1. 


(5.5) 


Combining  (5.1)  with  (5.4)  yields 


s*v  ds  +  (-fN)*V  ds  -  d 

~  dt  J 

y  y  y 


W(F)  dA, 


(5.6) 


suggesting  that  +fN^  be  viewed  as  a  fictitious  nominal  "shock  driving  trac¬ 
tion"  exerted  on  the  shock  by  the  surrounding  material;  the  scalar  f  deter¬ 
mines  the  magnitude  of  this  traction.  This  concept  is  related  to  the  notion 
of  the  "force  on  a  defect"  introduced  by  Eshelby  [9],  The  expressions 
(5.5),  (5.3)  together  comprise  a  special  case  of  a  result  given  in  [14]; 
see  also  Rice  [15], 


When  treated  from  a  thermomechanical  point  of  view,  the  dissipation 
rate  d(t)  can  be  shown  to  be  identical  to  the  product  of  temperature  and 
entropy  production  rate,  provided  that  the  temperature  is  spatially  uniform 
and  constant  in  time.  The  Claus ius-Duhem  inequality  then  requires  the  dis¬ 
sipation  rate  d(t)  to  be  non-negative  for  all  sub-regions  II  and  all 
instants  t;  see  [8].  Equivalently,  one  must  have 


f  Vn  >  0,  xect,  to  <  t  <  t]_, 


(5.7) 


where  Vn  is  the  normal  velocity  of  a  point  on  the  shock  Ct 


vn  “  V  •  N. 


(5.8) 


The  dissipation  inequality  (5.7)  is  trivially  satisfied  in  the  special 
case  when  f  happens  to  be  zero.  In  general,  however,  this  is  not  the  case, 
and  given  an  equilibrium  state,  (5.7)  specifies  the  direction  in  which  the 
the  shock  may  move  in  a  quasi-static  motion  commencing  from  this  state. 


WWW! 


Following  Yatomi  and  Nishimura  [13],  one  can  eliminate  the  shock 


orientation  N,  from  (5.5),  (5.3)  as  follows:  From  the  discussion  following 


(3.13), 


;F]]+  -  [[FN]]+®  N. 


(5.9) 


Using  this  and  (2.19)^  in  (5.5),  (5.3)  yields  the  desired  result: 


f  “  [[W]]+  -[[F]]+-  a. 


(5.10) 


where  the  dot  product  of  two  tensors  A  and  B  is  defined  to  be  A*B  « 

‘  A*  /W  S*m 


trace (AB^) . 


A  more  illuminating  expression  for  f  may  be  obtained  as  follows.  First, 


from  (5.5),  (5.3),  (2.19)  and  (3.13),  one  finds 


f  -  [ [W] ]  +  K  {a N  •  i)  (F~n  •  N) . 


(5.11) 


Using  (3.10),  (2.6)  in  (5.11)  provides 


f  -  [ [W] ]  +  (rn  •  H)k ; 


(5.12) 


rn*i  is  the  shear  traction  parallel  to  the  shock,  while  k  is  the  amount  of 


the  additional  shear  introduced  in  the  representation  (3.13) 


Finally,  from  (5.12)  we  derive  an  expression  for  f  that  is  appropriate 


for  plane  deformations  of  isotropic,  incompressible,  elastic  materials  and 


that  does  not  involve  the  shock  orientation.  By  (2.11)2,  (2.14)  and  (3.3) 


we  may  write  (5.12)  as 


WWW 
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*•1 


't'l 


m 


i 


a 


f  -  [ tW] ]  +  k  M(k+)  i-Gn. 


(5.13) 


i 

Using  (3.18)  to  express  k  in  terms  of  i«Gn  and  then  eliminating  the  latter 
from  (5.13)  with  the  help  of  (4.4),  we  find 


f  -  [[W]]+  +  (I  -  g2  -g-2)(M(k+)/M(k.))[ [M(k)]]+. 


(5.14) 


T  9 

Writing  I  -  2  +  k+z  and  using  (4.8)  to  eliminate  g  from  (5.14)  yields 


f  -  [ [W(2+k2) ] ] 


Equivalently,  we  may  write 


M(k+)  M(k.) 

-  i 

M(k+)  +  M(k.) 


(k+2  -  k.2). 


(5.15) 


f  -  H(k+1k_)  on  Ct  ,  tj  <  t  <  tj  , 


where  H  is  the  function  defined  on  E  by 


(5.16) 


H(k1,k2)  -  r (k)  dk  - 


r(k1)r(k2) 
r(k^)k2  +  r(k2)k^ 


(kx2  -  k22),  (k1,k2)eS. 


(5.17) 


By  (5.16),  the  shock  traction  f  at  a  point  on  the  shock  Ct  depends  only  on 
the  local  amounts  of  shear  k+,  k.  associated  with  the  deformations  on  the 
two  sides  of  Ct  ;  in  particular,  f  does  not  depend  on  the  the  orientation 
of  the  shock. 


Although  derived  here  for  plane  deformations,  the  expression  (5.16), 
(5.17)  for  the  shock  traction  is  formally  identical  to  the  corresponding 
result  derived  by  Yatomi  and  Nishiraura  for  anti-plane  deformations  (see 
equation  (3.3)  of  [13]).  This  is  not  surprising,  in  view  of  the  represen¬ 
tations  (2.3)  for  plane,  area-preserving  deformations. 


S'? 


--23-- 


A  A 

In  the  special  case  of  a  normal  shock,  we  have  r(k+)  -  r(k.)  by  (4.15), 
so  that  (5.16),  (5.17)  reduce  to 


f-  r  ( k )  dk  -  r  ( k+ )  (k+  -  k.). 


(5.18) 


From  (5.18),  the  value  of  f  in  the  case  of  a  normal  shock  mav  be  inter¬ 
preted  geometrically  as  the  difference  -  A2 ,  where  A-^  is  the  area  under 
the  stress-strain  curve  between  k_  and  k+,  and  A2  is  the  area  of  the 
rectangle  on  the  same  base  with  height  r(k+) . 


In  general,  H(k^,k2)-0  defines  a  curve  lying  in  the  region  2  of  the 
(ki , k2) -plane .  The  shock  traction  corresponding  to  points  on  this  curve  is 
zero,  and  therefore  so  is  the  local  rate  of  dissipation.  If  at  some 
instant,  the  limiting  deformation  gradient  tensors  F  and  F  at  a  particle  on 
Ct  are  distinct  but  such  that  H(k+,k_)  -  0,  we  say  that  this  particle 


is  in  a  "I 


i" .  For  a  particle  which  happens  to  involve  a  normal 


shock,  the  Maxwell  states  as  characterized  through  (5.18)  with  f  -  0  sat¬ 
isfy  the  well-known  "equal-area  rule". 


Figure  5  shows  the  curve  H(k^,k2)“0  for  a  material  characterized  by  a 
shear  stress  response  function  r(k)  of  the  form  shown  in  Figure  2.  For 
the  purpose  of  drawing  the  figure,  we  have  used  the  particular  response 
function  given  in  (4.14).  The  shaded  area  indicates  points  at  which  H  is 
positive;  the  unshaded  area  corresponds  to  negative  values  of  H.  When  *(k) 
has  the  form  shown  in  Figure  2,  one  can  show  that  the  maximum  and  minimum 
values  of  H  occur  on  the  boundary  of  the  set  2,  i.e.  at  certain  normal 
shocks . 
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